home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
IRIX 6.2 Development Libraries
/
SGI IRIX 6.2 Development Libraries.iso
/
dist
/
complib.idb
/
usr
/
share
/
catman
/
p_man
/
cat3
/
complib
/
BALANC.z
/
BALANC
Wrap
Text File
|
1996-03-14
|
3KB
|
67 lines
____BBBBAAAALLLLAAAANNNNCCCC((((3333FFFF)))) ____BBBBAAAALLLLAAAANNNNCCCC((((3333FFFF))))
NNNNAAAAMMMMEEEE
BALANC, SBALANC - EISPACK routine. This subroutine balances a REAL
matrix and isolates eigenvalues whenever possible.
SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee bbbbaaaallllaaaannnncccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh
ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaa((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssbbbbaaaallllaaaannnncccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, aaaa,,,, lllloooowwww,,,, iiiigggghhhh
rrrreeeeaaaallll aaaa((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
On INPUT
NNNNMMMM must be set to the row dimension of two-dimensional array parameters
as declared in the calling program dimension statement.
NNNN is the order of the matrix.
AAAA contains the input matrix to be balanced. On OUTPUT
AAAA contains the balanced matrix.
LLLLOOOOWWWW and IGH are two integers such that A(I,J) is equal to zero if
(1) I is greater than J and
(2) J=1,...,LOW-1 or I=IGH+1,...,N.
SSSSCCCCAAAALLLLEEEE contains information determining the
permutations and scaling factors used. Suppose that the principal
submatrix in rows LOW through IGH has been balanced, that P(J) denotes
the index interchanged with J during the permutation step, and that the
elements of the diagonal matrix used are denoted by D(I,J). Then
SSSSCCCCAAAALLLLEEEE (J) = P(J), for J = 1,...,LOW-1
= D(J,J), J = LOW,...,IGH
= P(J) J = IGH+1,...,N. The order in which the interchanges
are made is N to IGH+1, then 1 TO LOW-1. Note that 1 is returned for IGH
if IGH is zero formally. The ALGOL procedure EXC contained in BALANCE
appears in
BBBBAAAALLLLAAAANNNNCCCC in line. (Note that the ALGOL roles of identifiers K,L have been
reversed.) Questions and comments should be directed to B. S. Garbow,
Applied Mathematics Division, ARGONNE NATIONAL LABORATORY
PPPPaaaaggggeeee 1111